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Abstract: Inverse nodal problem on diffusion operator is the problem of finding the potential 
D ' functions and parameters in the boundary conditions by using nodal data. In particular, we solve the 
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reconstruction and stability problems using nodal set of eigenfunctions. Furthermore, we show that the 
space of all potential functions q is homeomorphic to the partition set of all asymptotically equivalent 
nodal sequences induced by an equivalence relation. To show this stability which is known Lipschitz 
stability, we have to construct two metric spaces and a map &dif between these spaces. We find 
that Qdif is a homeomorphism when the corresponding metrics are magnified by the derivatives of q. 
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Basically, this method is similar to pQ and [2] which is given for Sturm-Liouville and Hill operators, 
respectively and depends on the explicit asymptotic expansions of nodal points and nodal lengths. 
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CN . 1. Introduction 

Inverse spectral problems is about recovering operators by using their spectral charachteristics 
as spectrum, norming constants and nodal points. Such problems have been an important research 
issue in mathematics and have many applications in natural sciences. Inverse spectral problems 
are divided into two parts. One of them is inverse eigenvalue problem and the other one is inverse 
nodal problem. Inverse eigenvalue problems have been studied for along time by many authors 



EL ED) EL E], EL E]. Inverse nodal problem was first posed and solved by McLaughlin [9] who 
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showed that the knowledge of a dense set of nodal points of the eigenfunctions alone can determine 
the potential function of the Sturm-Liouville problem up to a constant. Independently, Shen and 
his coworkers studied the relation between the nodal points and the density function of the string 
equation [10J. Recently, many authors have studied inverse nodal problem for different operators 

m, mi, mi, mi, mi, mi- 

Consider the following boundary value problem generated by the differential equation 

Ly = y" + [A 2 - 2Xp(x) - q(x)] y = 0, x G (0, vr), (1.1) 
and boundary conditions 

j/(fl)-hy(fi) = 0, (1.2) 

y'(7T)+Hy(7T) = 0. (1.3) 

where A is spectral parameter and h,H £ R; q £ W\ [0,tt],jj E Wf [0,7r] [T7j, [TB]. Since L 
is determined by its potential functions, we identify L with L pq . Problem (1.1)-(1.3) are called 

nnarlratir 1 npnnilt: nf ^rhrnrlinow nrM^rainr (nr riifFiisinn nripratnr^ S-lnmp vpr^inns; nf piffpnvaliip 



Let A n be the n— th eigenvalue, y(x, A n ) the eigenfunction corresponding to the eigenvalue A n 
and < < x^ < ... < < tt be the nodal points of the re— th eigenfunction y{x, A n ). In 
otherwords, A n ) = 0, j = 1, 2, n — 1. Also let Ij = [xj ,xj^ x ] be the j— th nodal domain 

of the n— th eigenvalue and let 1^ = x < ^ 1 — be the nodal length. Define Xq = 0, x™ = it. We 
also define the function j n (x) to be the largest index j such that < x™ < x. Thus, j = j n (x) if 
and only if x € [x™,x™ +1 ). Define X = {x™},re > 0, j = 1,2, ...n — 1. X is called the set of nodal 
points of (1.1)-(1.3). 

Inverse nodal problem for diffusion operator is to determine potential functions and parameters 
in the boundary conditions. This type problems have been studied by many authors [27], [55], [2S], 

[an], m- 

This study is organized as follows: in section 2, we mention some physical and spectral properties 
of diffusion operator. In section 3, we give a reconstruction formula for potential function and some 
important results for the problem (1.1)-(1.3). Finally, we solve Lipschitz and high order Lipschitz 
stability problems for diffusion operator in sections 4 and 5, respectively. 

2. Some Physical and Spectral Properties of Diffusion Equation 

Jaulent and Jean [32] stated the actual background of diffusion operators and discussed the 
inverse problem for the diffusion equation. Also, Gasymov and Guseinov studied the spectral 
theory of diffusion operator [T7] . 

The problem of describing the interactions between colliding particles is of fundemental inter- 
esting in physics. In many cases, a description can be carried out through a well known theoretical 
model. In particular, one is interested in collisions of two spinless particles, and it is supposed 
that the s— wave binding energies and s— wave scattering matrix are exactly known from collision 
experiments, s— wave Schrodinger equation with a radial static potential V can be written as 

y" + [E — V(x)\ y = 0, x>0 (2.1) 

where the potential function depends on energy in some way and has the following form of energy 
dependence 

V(x,E) = U(x) + 2VeQ(x). (2.2) 

U(x) and Q(x) are complex-valued functions. (2.1) reduces to the Klein-Gordon s— wave equation 
with the static potential Q(x), for a particle of zero mass and the energy yE with the additional 
condition U(x) = -Q 2 (x) [32]. 

The Klein Gordon equation is considered one of the most important mathematical models in 
quantum field theory. The equation appears in relativistic physics and is used to describe dispersive 
wave phenomena in general. In addition, it also appears in nonlinear optics and plasma physics. 
The Klein- Gordon equation arise in physics in linear and nonlinear forms [33 J. 

Now, we will consider Klein Gordon wave equation. After some straightforward computations, 
it turns diffusion equation which is given in (1.1). This form of Klein Gordon equation was first 
given in [33]. But here we shall improve our understanding. Let consider following Klein Gordon 
wave equation 

r / Ft \ 2 / 1 -x\ 2 ~\ 



This equation represents a spinless particle of charge e and mass m in a scalar potential (f> and 
vector potential A (r, t) where natural units h = c = 1 have been used. In case of ~t(r,t) = and 
the scalar potential to be time independent, equation (2.3) is reduced the following form 



(2.4) 



If we set \l/(r, t) = ip(r)e tEt , equation (2.4) takes the form 



where 



V 2 + (E-e(t>) 2 <p(r) = m 2 <p(r), (2.5) 

tf t (r,t) = -ip(r)iEe- iEt 
V tt (r,t) = -^(r)£ 2 e-^. 

This equation is a well-known form of Klein-Gordon equation. The wave equation (2.5) with a 
spherically symmetric potential energy may be written in spherical coordinates 



r 2 dr V dr J r 2 sin 9 89 



sine— + 
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89 J ' r 2 sin 2 i 



+ (E - ecf>) } f = m 2 if 



(2.6) 



We first separate the radial and angular parts by substituting 



<p{ r ,e,<f>) = R(r)Y{e,, 

into equation (2.6) and dividing through by tp, we get 



LA. 



,dR\ 
dr J 



(E - ecj)) 2 - m 2 



1 

Y 



i a 

sin 9 89 



sin ( 



dY 
~89 



+ 



1 8 2 Y 
sin 2 9 ~8^_ 



(2.7) 



Since the left side of equation (2.7) depends only on r, and the right side depends only on 9 and 
<fi, both sides must be equal to a constant that we call 1(1 + 1) where I is orbital quantum number. 
Then, the equation (2.7) gives us a radial equation with the scalar potential V = ecfr for the Klein 
Gordon equation 
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,dR\ 



r 2 dr V dr ) 



where I = 0, 1, 2, ... Substituting i? = 



(£ 2 - 
V>(r) 



m 



! ) +V(V- 2E) 



1(1 + 1) 



R = 



in equation (2.8), we obtain 



(2.8) 



<(r) + 



iT 2 + V(V - 2£7) - 



/(/ + !) 



V>(r) = 0, 



where -E 2 = K 2 + m 2 . In case of / = m = 0, we get 



</>"(r) + \V 2 - 2KV] Mr) = -K 2 ib(r). 



(2.9) 



There are many important spectral properties of eigenvalues and eigenfunctions for the problem 

(1.1) -(1.3). We collect some of these in the lemmas below. 

Lemma 2.1. The eigenvalues of the operator L pq are simple. 
Proof: It can be easily proved by using similar way with |35j . 

Lemma 2.2. The operator L p ^ q defined by (l.l)-(l.S) is in fact symmetric on space L2[0, 7r]. 
Proof: Let u and v be twice differentiable functions which satisfy the boundary conditions 

(1.2) -(1.3). This lemma can be esily proved by using integration by parts on these two functions u 
and v. Considering the definition of inner product on L2[0,7r], we get 

TT TT 

< L Pj g [u] , v >= J L ps [u] vdx = J {—u" + [q(x) + 2Xp(x)] u] vdx 



TT 

= hu(0)v{0) + v(tt)Hu(it) - v(tt)Hu(tt) - hu{0)v(0) + J {-v" + [q(x) + 2Xp(x)] v} udx 

o 

= < U, L P:q [v] > . 



Lemma 2.3. The eigenfunctions corresponding to different eigenvalues are orthogonal for the 
operator L Ptq . 

Proof: Let 4> m and <f> n be the eigenfunctions corresponding to the eigenvalues X m and A n , (X m ^ X r 
respectively. If <f) k denotes the k— th eigenfunction, we can integrate the identity 



to obtain 



-KA + $'j>m + 2 P(, x )4>m<Pn ( A m ~ K) = (A£, - A^) <p m 4> n , 



\2p(x) - X m - X n ] 4> m <j) n dx = 0, 



from which the result follows. 



Lemma 2.4. [TT], [30] Let X n ,n £ Z - {0} be the spectrum of the problem (l.l)-(l.S). It is 
well known that the sequence {X n : n = 1, 2, ...} satisfies the following asymptotic formula 



X n = n + c -\ 1 '—, 

n n 



where 



If 1 

co = - / p(x)dx,a = - 

TT J 7T 




h + H + 



1 



[q(x) + p 2 (x)] dx 



, |ci jn | 2 < oo. 



n=l 



Let consider the equation (1.1) with the initial conditions 

2/(0) = 0, y'(0) = 1, 



(2.10) 



(2.11) 
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Lemma 2.5. [27] The solutions of the problems (l.l)-(l.S) and (1.1), (1.3), (2.11) have the 
following forms, 

X 

h f sin \X (x — t)] 
(p(x, A) = cos (Xx) - - sin (Ax) + / [q(t) + 2\p(t)] <p(t, X)dt, (2.12) 

o 

x 

sin (Ax) f sin [A (x — t)} r , . „ , , ., , , , . , /„..„.. 
*P(x, A) = ^ + J >-± [q(t) + 2Xp(t)) V(t, X)dt, (2.13) 

o 

Lemma 2.6. [27] Suppose that q S Li[0,7r] and p G Wj^O, 7r]. Then, the nodal points and 
nodal lengths has the following asymptotic forms, as n —> oo 
(a) For f/ie problem (l.l)-(l.S), 



x] 



J~2 I* 



3 

^ + ^2 y [1 + cos (2A n i)] [g(t) + 2A n p(t)] dt + o (J^j , (2.14) 



^ = ^ + 7772 / [1 + cos (2A n t)] [g(i) + 2A n p(t)] + o f ^ V ( 2 - 15 ) 
A„ 2A n y \ A n/ 

(b) For the problem (1.1), (1.3), (2.11), 

X 3 

^ J [1 - cos (2A n t)] [g(t) + 2A n p(f)] <ft + o (J-^J , (2.16) 



n _ jlT 1 

Xj " a; + 2~x* 



o 

n 



/ i=T^ + ^ / [l-cos(2A„t)][?(i) + 2A n p(t)]dt + o(i j. (2.17) 

Lemma 2.7. Suppose that q G Li[0, 7r]. Then, for almost every x G (0,7r) , u>ii/i j = j n [x) 

lim A n / q(t)dt = q(x). 



n— >oo 
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Proof: It can be proved by using similar way with pQ. 

The above equalities are still valid even if x™ and x™ +1 are replaced by Xj and Xj +1 . 

3. A Reconstruction Formula for Potential Function and some Important 
Results 



In this section, we will give a reconstruction formula for the potential function of diffusion 

nnpratnr nnrl snmp rpsnilt« TVipn wp «Vinll crivp rlpfinitinns nf rl^ nnrl dr. A^/p will rlpnntp T J} ( crr\A 



the space Qdif as a collection of all diffusion operators L VA and the space S^/ as a collection of 
all admissible double sequences of nodes such that corresponding functions are convergent in L±. 
A pseudometric ds dl/ will be defined on Essentially, d-£ dif (X,X) is so close to 



n-l " 

d (X, X) = lim n 2 7r ^ | L£ - TT k | + n / \p(x) - p(x) \ dx, 

k=l n 



(3.1) 



where ^dif, L' k = X k+1 — X% and h\ is defined similarly. 

If we define X ~ X if and only if ds dif (X, X) = 0, then ~ is an equivalence relation on 
T<dif and ds di/ would be a metric for the partition set £^ = E^// ~ . Let S^i/i C E^i/ be the 
subspace of all asymptotically equivalent nodal sequences and let = ^<fi/i/ ~ • Let <I>^ y be 

the homeomorphism between the spaces O^/ and We call $dif as a nodal map for diffusion 

operator. 

Theorem 3.1. Let q G Li[0, 7r] and p G W 2 2 [0, 7r]. Define F n by 



F„ = n 



n-l 



^ 2n 2 L£ - 2nvr - 2p(rc) 



.fc=i 



for the problem (l.l)-(l.S). Then, F n converges to q poinwisely almost everywhere and also in L\ 
sense. Furthermore, pointwise convergence holds for all the continuity points of q. 

Proof: We shall consider the asymptotic formulas for the nodal lengths of the problem (1.1)- 
(1.3). Observe that by Lemma 2.6, we have 



TT 



1 



"i+1 



2X n 7T 



[1 + cos (2X n t)\ [q(t) + 2X n p(t)} dt + o 



1 



"i+i 



And, after some algebraic computations, we get 

2A 2 (l]X n - tt) - 2X nP (x) = X n J q(t)dt + o(l). 

If we consider the left side of this equality with the asymptotic formula of A n , it yields 
2A 2 (qX n - tt) - 2X nP (x) = X n {2X 2 n q - 2X n n - 2p{x)) 



n + o | — 

n 



n + o ( — 

n 



q-2 



n + o ( — 

n 



TT — 2p(x) 



= n (2n 2 q - 2nir - 2p(x)) + o(l). 
Hence, to prove Theorem 3.1., it suffices to show Theorem 3.2. (b). 
Remark: F n can be obtained similarly for the problem (1.1), (1.3), (2. 11). 
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(a) For the problem (l.l)-(l.S), 



h = lim 2X n n ( j - - - -^XJ ) . 



1 _ An 

2 vr 



(b) For almost every x E [0, ir] 

q(x) = lim 2A n [A^/™ - A n 7r - p(x)] . 

where q £ Li[0, 7r]. 
Proof: 

(a) From the Lemma 2.6, we know the following asymptotic formula for the problem (1.1)-(1.3), 

(j - - V x " 

X J = \J + w I [1 + cos {2xnt)] m + 2xnp{t)] dt + ° (f) • (3,2) 

After some computations in (3.2), we get 

x i 

Xn7T ^ X i) = ^~lj [1 + COS {2Xnt)] [q{t) + 2XnP{t)] dt + ° (i) • 

o 

Since X™ goes to zero as n — > oo, we obtain 

1 A 



h = lim 2A n vr I j - - - ^X? 
(b) Part (b) can be proved by using similar procedure with in [27] . 

Theorem 3.3. Given L Pj(? in fldif, the set of nodal points {x^} is also asymptotically nodal 
to L Pt g itself. 

Proof: We can easily prove by using similar way with in pQ. 

4. Lipschitz Stability of Inverse Nodal Problem for Diffusion Operator 

In this section, we solve a Lipschitz stability problem for diffusion operator. Lipschitz stability 
is about a continuity between two metric spaces. To show this continuity, we will use a homeomor- 
phism between these two spaces. Stability problems were studied by many authors p], [2], |36j, 
[37], [38]. 

Definition 4.1. Let N' = N — {1, 2}. 

(i) 

^dif = {q £ £i(0, vr) : q is the potential function of the diffusion equation} 
T>dif = The collection of the all double sequences defined as 



for each n £ N. 

(ii) Let X e and define X = {X%} where L™ = X™ +1 - X£ and I£ = (X™, X% +1 ) . We say 
X is quasinodal to some q G f^i j if X is an admissible sequence of nodes and satisfies (I) and (II) 
below : 

(I) X has the following asymptotics uniformly for k, as n — > oo 



X n = \ U_ + l — ),k = l,2,...,n 

n \n z 



A fc — " 

and the sequence 



n-l 



^ 2n 2 L£ - 2nvr - 2p(a 



fc=l 



converges to q in Li for the problem (1.1)-(1.3). 

(II) X has the following asymptotics uniformly for k, as n — > oo 

= — + o(^Yfc = l,2,...,n. 
n \n z J 

Similarly, for this case F n converges to q in L\ for the problem (1.1), (1.3), (2.11). 

Definition 4.2. Suppose that X, X € E^/ with L£ and as their respective grid lengths. 

Let 

n-l * 
S„(X,X) =n 2 ^^|L^-Lfc|+n J \p(x) - p(x)\ dx . (4.1) 



fe=l 

Define 



s n (x,x) 



d {X,X) = lim S n (X,X) and cfe (X,X) = lim - n v ' ' . 
This definition was first made by PQ for Sturm-Liouvile operator. Since the function /(x) 



is monotonic, we get 



1 + x 

. _ d (X,X) r , 
X,X) = — v / L= G 0,vr . 
dl/ 1 ' ; l + d (X,X) 1 ' J 

Conversely, 



d (x,x)= d *«f( x ; x ) 



This equalities can be obtained easily. 

Lemma 4.1. Let X, X e 

(a) ds dif is a pseudometric on S^j. 

(b) 7/ X and X belong to different cases, then d^ diJ (X, X) = 1. 

(c) // X belongs to case (I) or case (II), then 

L k = - + 0[^-] ,k = l,2,...,n. 



i) The interval 5 n ^ between the points XJ? and X\ has length 0(— ), as n — > oo. 

ii) For all x £ (0, 7r), define J n (x) = max{£; : -X'i™'' < x} so £/ia£ A; = J n (x) if and only if 
x £ [Xj,Xj +1 ). Then, | J n {z) — J n {z)\ < 1 for sufficiently large n. 

Proof: It can be proved similar to [1], [2]. 

After following theorem, we can say that inverse nodal problem for diffusion operator is Lipschitz 
stable. 

Theorem 4.1. The metric spaces (Qdif, H-Hi) an d (T^difi/ ~,c?s di /) o- re homeomorphic to each 
other. Here ~ is the equivalence relation induced by d-£ dif . Furthermore 

2d Sdif {X,X^ 
llq qh l-dx dif {X,X)' 

where d^ dif (X,X) < 1. 

Proof: In view of Lemma 4.1., we only need to consider when X, X £ ^dif belong to same 
case. Without loss of generality, let X, X belong to case I. We have to show 

\\q-q\\ 1 = 2d {X,X). 

According to the Theorem 3.1., F n and F n convergence to q and q, respectively. If we use 
definition of norm on L\ for the potential functions, we get 



?lli < 2n 3 



dx + 2n J \p{x) — p(x)\ dx + o(l) 
n o 



and after some algebraic operations 



\q - q\\i < 2re 



t n j n 

L J n {x) ~ L Jn{x) 



dx + 2n" 



L J n (x) ~ L J n (x) 



dx + 2n j \p(x) — p(x)\ dx + o(l). 

(4.2) 



Here, the integrals in the second and first terms can be written as 



-yn -yn 

L J n (x) ~ L J n (x) 



and 



n-l 



L J n {x) ~ ^Jn(x) 



dx ~ ~ ^ \ L k ~ L k \ i 



k=l 



respectively. If we consider these equalities in (4.2), we get 



1 1 II 



f 



and 



n-1 

q-q\\ 1 <2n 2 TrJ2\ L k- Tl k\+2n J \p(x) - p(x)\ dx + o{\). (4.3) 



k=l 

Similarly, we can easily obtain 

n-1 



q-q\\ 1 >2n\J2\ L k- 1 %\ +2n J \p{x) -p{x)\dx + o{\). (4.4) 



k=l 

Considering (4.3) and (4.4) together, it yields 



q-q\\i =2n 2 TT^\Ll-L n k \+2n / |p(x) - p(x)\ dx. 



o 



fc=i 

The proof is complete after by taking limit as n — > oo. 
5. High order Lipschitz stability 

In this section, we will solve a high order stability problem for diffusion operator. For n £ N, 

^dif(N) = {q€L 1 [0,vr] : q G C* +1 [0,vr]} . 

It has been proved in [31] that m— th derivative of potential function q for Diffusion operator 
can be approximated by some difference quotient of nodal length 5 m Z™ where 5 m is m— th order 
difference quotient operator defined as [T] 

fa "> = !2±i_^ and jm ») = 5 a i , 

J ml J in) 

Note that 5 and 5 m operators depend on the double sequence {a[ n ^}. 

Let T,dif(N) be the set of asymptotically equivalent nodal sequences in Qdif(N). Let dfi diJ (jv) 
and -D£ dl/ (jv) be some metrics on fi^(JV) and E^y(iV) respectively magnified by the Li norms of 
derivatives of the potential functions. We find that the nodal map $>dif is still a homeomorphism 
under these strengthed metrics. 

If we define X ~jv X if and only if D^, dif (N) [X, X) = 0, then ~tv is an equivalence relation on 
^dif(N). Hence dn di/ (jv) is a metric on 0, dif (N). 

Definition 5.1. Let I,Ie Y, di f and X£ = X£ = 7r for > n. For m = 1, 2, TV, let 

1 n-m-2 n-m-l 

S m , n (X,X) = XI l^Ll-^l+Xn I \6 m p(x 7 )-5 m p(xj)\dx 

+A„ / p^ m \x) -p (m \x) dx. (5.1) 



Define 



d m (X, X) = lim S m „ [X, X) and dv^rmi (X.X) = lim '■ — ==-. 

We have to use the following reconstruction formula for the function in Theorem 5.1. to 
prove Theorem 5.3. 

Theorem 5.1. |31j Let q is real valued, (n + l)—th order continuous function from the class 
Li[0,7r] for N > 1 in (1.1), and let j = j n (x) for each x G [0, ir]. Then, as n — )• oo 



q(x) = X n [2X 2 J] - 2vrA n - 2p(x)] + O (-) , 



and, for all m = 1,2, N 



2A 



3/2 



q M {x) = _ 2X n 5 m p(x J ) - 2X nP ( m \x) + O(l). 

7T 

Theorem 5.2. The metric spaces (tt dif (N),d nMf(N) ) and (S difl (N)/ ~N,Ds dif (N)) are 
homeomorphic to each other, where ~_/v is the equivalence relation induced by D^ dif ^y 

Proof: In view of Lemma 4.1. we only need to consider X, X £ ^dif belong to same case. 
We shall show that \\q^ -q^\\ = 2d m (X,X). Hence, to prove Theorem 5.2., it suffices to prove 
Theorem 5.3. 

After Theorem 5.3., we can say that the inverse nodal problem for the diffusion operator is high 
order Lipschitz stable. 

Theorem 5.3. Suppose that q and q are both Nth order continuous functions from the class 
L\[0, 7r]; L Pt q,Lp^g belong to same case. Let X and X be the corresponding asymptotically equivalent 
nodal sequences. Then, for all m = 1,2, ...,N 

q (m)_q(m)j =2d m (X,X). 
Proof: By Theorem 5.1 and Lemma 2.6, we can write 



2A 



3/2 



q (m) _-(m) = pn L n _ gnq*} + ^ [fi^fa) ~ S m P (xj)]+2X n (x) - (x)J + O (1) , 

and 

2 \3/2 _ 

, ,,, ,„ .-,„r'n r ,-,„-. -i . pM( x ) +0 (1) 



_ g( m ) 



7T 



Then, by using the definition of norm on L\ space, we obtain 



_ 7?( m ) 



I U m \x) -q( m \x)\dx, 



and 





n — m — 1 


q( m ) — q( m ) 













2\T 



IT 



[ S m L n _ 6 m L n^ + ^ [ § ™p( Xj ) - 5 m p(xj)] 



+2A r 



p< m \x) -p {m \x)] +o(l) 



(5.2) 



when n is sufficiently large. By using the property of triangle inequality in (5.2), we get 

i 

\S m p(x 7 ) -5 m p(xj)\dx 



9X 3/2 



< 

1 7T 



+2A T 



LTL} - 5 m L7 dx + 2A„ 



Xf 



-V': 



p^ m \x) - p {m \x) dx + o{l), 



Xi 



and after some computations, 



? M _ q(™) 



< ?*= 

1 7T 



3/2 



| (5 m L n_ (5 m L «| dx + 



2 A 



3/2 



7T 



\S m Lj - 5 m Lj\ dx 



X" 



xi: 



+2A n J 5 m p{x 1 )-5 m p{x J )dx + 2\ n j \p im \x)-p i:m) {x)\dx + o{l) 



Xi' 



X',' 



Now, by using [[39 , Lemma 2.1] and Lemma 4.1., we obtain 



x '■ 



n—m—2 



J \5 m L n j-5 m L n j\dx= \5 m Ll-5 m L n k \\Ll\, 



and 



Vfl- 

jY n-m-l 



\S m L} - 5 m Lj\ dx = 



n' 



x™ 



Moreover, L™ = ^ + O . Thus, 



2A 



3/2 



X'; 



7T 



\S m L n j-5 m Lj\dx < 



2\T 



7T 



X! 



1 x ox 3/2 n-m-2 

re 4 / 7r 

7 fc=l 



x n-m-2 

J_j + 2n 1 / 2 £ |^-<5 m I"| + ^n 3 / 2 nO(n- 2 ) 
fc=i 



Y \5 m L n k -5 m L n k \ 



nn \ 1 1 - 



Therefore, 



n—m—2 



< 2nV 2 \S m Ll-5 m L n k \+2X n J \5 m p(x 7 ) - S m p(xj)\ dx 

xv- 



k=i 



n — m—1 

+2A n J Y m \x)-p i - rn \x)\dx + o(l). 



(5.3) 



If we take limit as n — > oo and use the metric definition, we obtain 

< 2d m (X,X). 



q( m ) _ g( m ) 



Using similar procedure, 



q( m ) _ q( m ) 



(5.4) 
(5.5) 



>2d m (X,X). 
Then, considering (5.4) and (5.5) together, we get 

| g (m)_ 5 (m)|| =2d m (X,X). 

This completes the proof. 
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